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The relationship between magnetic reconnection and plasma turbulence is investigated using mul¬ 
tipoint in-situ measurements from the Cluster spacecraft within a high-speed reconnection jet in the 
terrestrial magnetotail. We show explicitly that work done by electromagnetic fields on the particles, 

J -E, has a non-Gaussian distribution and is concentrated in regions of high electric current density. 

Hence, magnetic energy is converted to kinetic energy in an intermittent manner. Furthermore, we 
find the higher-order statistics of magnetic field fluctuations generated by reconnection are char¬ 
acterized by multifractal scaling on magnetofluid scales and non-Gaussian global scale invariance 
on kinetic scales. These observations suggest J ■ E within the reconnection jet has an analogue in 
fluid-like turbulence theory in that it proceeds via coherent structures generated by an intermittent 
cascade. This supports the hypothesis that turbulent dissipation is highly nonuniform, and thus 
these results could have far reaching implications for space and astrophysical plasmas. 


Introduction .—Turbulence is a universal nonlinear phe¬ 
nomenon that is ubiquitous in space plasmas [ij. It pro¬ 
duces a cascade of coherent structures in neutral fluids 
0 and plasmas d, Q . These are concentrated structures 
that are phase correlated over their spatial extent and 
have relatively long lifetimes, such as current or vorticity 
sheets. Indeed, current sheets have been observed ex¬ 
tensively in turbulent plasmas, and are associated with 
magnetic reconnection in the solar wind [5|, at the mag¬ 
netopause [ 1 ] , and in the magnetosheath [3, 0 . In the 
terrestrial plasma sheet, in-situ coherent structures dis¬ 
play signatures of intermittent turbulence in the 

form of rare large amplitude fluctuations that are highly 
non-Gaussian. These spatially inhomogeneous turbulent 
flows have been proposed as central to plasma sheet dy¬ 
namics [IMl, and thus critical to understanding how 
stored electromagnetic energy in the magnetotail is con¬ 
verted into plasma energy. Here we consider whether 
coherent structures generated by magnetic reconnection 
reflect the nonlinear dynamics of intermittent turbulence 
and might be sites of nonuniform dissipation. These long¬ 
standing questions [1^ are the subject of this Letter. 


The nature of turbulent dissipation within collisionless 
plasmas remains an open problem. A strong turbulent 
cascade is far from equilibrium and smaller scale behavior 
is driven by larger scale dynamics, with faster response 
times for decreasing scales 15[ . In addition, intermittent 
turbulence generates small-scale coherent structures that 
are responsible for nonuniform dissipation (Sin. For 
neutral fluids, the Kolmogorov refined similarity hypoth¬ 
esis (hereafter KRSH; BEi) relates the statistics of 
increments of the velocity field on a given spatial scale to 
local averages of the dissipation rate on the same scale. 
Hence, large intermittent fluctuations on small scales are 
concomitant with high local concentrations of dissipa¬ 


tion. While KRSH is unproven, it is well supported 
in hydrodynamics and lies at the heart of modern 
fluid turbulence theory. In contrast, KRSH lacks verifica¬ 
tion and even precise formulation for collisionless plasmas 
since these introduce significant complications that have 
not yet been overcome. Amongst these complications is 
the inability to write an explicit form of the dissipation 
function, which is well known in viscous hydrodynam¬ 
ics and visco-resistive magnetohydrodynamics. However, 
plasma turbulence is well described by ideas that parallel 
its fluid antecedents, and thus it is instructive to test the 
hypothesis that coherent structures are linked to nonuni¬ 
form dissipation using, as a surrogate measure, the work 
done by the electromagnetic fields, J-E in an appropriate 
reference frame. 

Analysis .—We use 450 Hz burst mode magnetic field 
B and electric field E measurements from the FGM 2l| , 
STAFF 221 and EFW [2^ instruments onboard the Clus¬ 
ter spacecraft. While components of the dc electric field 
in the spacecraft spin plane are measured directly, the 
third component is reconstructed assuming E • B = 0. 
We also use 4 s resolution proton moments from the CIS 
CODIF experiment 241. 

The curlometer technique is used to estimate 

the current density J through a tetrahedron formed by 
four spacecraft, where the Maxwell-Ampere law is writ¬ 
ten as: 


Uo^ijk ‘ (Ar^/e X Ar^/i;) — ABj/^ ■ Ar^/^ ABj^, ■ Ar^/g (1) 

where i, j and k are the spacecraft indices, is the 
average current density normal to the surface made by 
spacecraft i, j and k, Ar^fc = is the distance 

between spacecraft i and k, and AB^fc = B^ — B^ is 
the magnetic field difference between spacecraft i and k. 
The total average current density is then determined by 
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projecting the current normal to three faces of the tetra¬ 
hedron into suitable cartesian coordinates. However, this 
technique is not without its limitations 27]. The main 
assumptions are that the spatial variation of the mag¬ 
netic field is a linear function of the spacecraft separation, 
such that J is constant over the tetrahedron, and that the 
medium is stationary. Since non-stationarity leads to the 
generation of nonlinear gradients, the only source of er¬ 
ror to consider is the nonlinear variation of the magnetic 
field. This is determined by computing V • B from: 

V • B|Arifc • Ar^-fe x Ar^zl = | ^ AB^^ ■ Ar^fc x (2) 


While B is solenoidal, the expression above can produce 
non-zero values which result from nonlinear gradients 
that are neglected in the estimate. Hence, V-B/ |V x B| 
provides an indicator of the error on curlometer estimates 
of J. We require this error to be less than 10%, and re¬ 
move all data that does not satisfy this condition. Note 
that certain spacecraft configurations and separations 
can reduce the accuracy of the curlometer technique, but 
these effects are likely insignificant in this study since the 
Cluster quartet are in a regular tetrahedral configuration. 

Results. —Figure 1 shows selected plasma and magnetic 
field data for a 12 min interval encompassing an earth¬ 
ward magnetic reconnection jet observed in-situ on 17 
August 2003 [e.g. 2^31|. The Cluster quartet are in 
the magnetotail lobe prior to 16:55 UT, and then enter 
the plasma sheet whilst located around [—16.8,5.6,3.2] 
Earth radii in geocentric solar magnetospheric (GSM) 
coordinates with spacecraft separations near 220 km. A 
high speed earthward flow (14 exceeds 1200 kms“^) is 
detected in the plasma sheet, where the number density 
n « 0.25 cm“^ and plasma beta /3 = nksT//2 ^,q) ss 
1 have typical values. However, the proton temperature 
is significantly elevated within the earthward flow, which 
is suggestive of proton heating by magnetic reconnection. 
The spacecraft exit the reconnection jet and enter the 
central plasma sheet at 17:03 UT. 

A tailward flow was detected by the spacecraft ahead 
of the earthward flow from 16:33 to 16:52 UT, and the 
associated was mostly negative (positive) for the tail- 
ward (earthward) flow. These correlated changes in 14 
and Bz are consistent with a tailward retreating X-line 
being swept past the spacecraft, and suggest that the 
earthward flow is close to this X-line. The negative sign 
of By PS —15 nT (roughly 60% of the asymptotic mag¬ 
netic field) agrees with the expected Hall magnetic field 
polarity in the southern hemisphere, eastward of the X- 
line. Indeed, the reversal in observed around 16:55:12 
UT could be associated with a reversal in the nearly field- 
aligned currents which close the Hall currents across the 
reconnection separatrices. Hence, the spacecraft may 
have observed an ion diffusion region with a moderate 
guide field. 



FIG. 1. An overview of a high speed earthward magnetic 
reconnection jet. The parameters from top to bottom are: 
GSM components of magnetic field and solar wind velocity, 
proton number density, proton temperature and current den¬ 
sity. Vertical dashed lines bracket the reconnection jet. The 
observations before and after the earthward jet are from the 
magnetotail lobe and central plasma sheet respectively. 


There are large fluctuations in the magnetic field and 
current density associated with the high speed reconnec¬ 
tion jet. Here we investigate the statistical properties of 
these fluctuations to identify signatures of intermittent 
turbulence. In order to determine the higher-order scal¬ 
ing of magnetic field fluctuations, the absolute moments 
of the increments 5B{t, r) = B{t-\-T) — B{t) are computed 
for each vector component B —>■ B^ , By or B^ ■ The mth 
order structure function is given by: 

1 ^ 

S^{r) = -Y.\^B{U^^)r (3) 

i=l 

where r is the time lag and N is the signal sample size. 
The higher-order structure functions progressively cap¬ 
ture the more intermittent fluctuations. These represent 
the spatial gradients responsible for dissipating energy in 
fluid-like turbulence. We examine the powerlaw scaling 
behavior of structure functions such that: 

S'^ir) oc (4) 

and C,{rn) are the scaling exponents. Figure 2 shows the 
GSM cc-component magnetic field structure functions and 
scaling exponents for the earthward flow data interval 
indicated in Fig. O- Note that the three magnetic field 
components exhibit essentially identical statistical scal¬ 
ing. The inertial and dissipation ranges are well defined 
with a break around 3 s, in agreement with the temporal 
signature of the proton gyroradius (~ 260 km). This sug¬ 
gests Taylor’s hypothesis [s^ is valid within the magnetic 
reconnection jet, and thus each spacecraft time series can 
be considered a spatial snapshot of the plasma. However, 
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this cannot be confirmed without knowing the character¬ 
istic timescale on which the observed fluctuations vary. 
Note that hereafter we assume time lags and frequency 
spectra are equivalent to spatial lags and wavenumber 
spectra. In effect we rely on some form of random sweep¬ 
ing of small scales [33 1. 

The scaling exponent C(2) is directly related to the 
power spectrum spectral index a when Eq. (4) is sat¬ 
isfied: C(2) = a — 1 [ 3 ^. This relationship only ap¬ 
plies to a < 3 when using two-point structure functions 
|4ll| . Here the inertial range follows a power law with 
a = 1.65 ± 0.03 which then steepens to a = 2.52 ± 0.02 
in the dissipation range. These are similar to other re¬ 
ported observations in magnetic reconnection diffusion 
regions 3^ 3^ and in the turbulent solar wind 371. 






FIG. 2. The magnetic field (a) structure functions and (b) 
inertial range scaling exponents. These structure functions of 
order 1-5 have been shifted along the vertical axis to facili¬ 
tate comparison of gradients. Linear fits for the inertial and 
dissipation ranges are also shown. The dissipation range (c) 
scaling exponents and (d) PDFs rescaled by their standard 
deviations. A Gaussian fit (dashed curve) is also applied. 


Figure 2 shows results from the higher-order scaling 
analysis of magnetic field fluctuations in the reconnection 
outflow. The errors on (^{m) are estimated as the sum of 
the regression error from Fig. 2(a) and the variation in 
C (m) found by repea ting the regression over a subinterval 
of the scaling range |38l |. For inertial range fluctuations, 
(m) is nonlinear in m, which is typical of hydrodynamic 


[161] and magnetofluid 17| turbulence. This behavior is 
associated with a statistical distribution of energy dis¬ 
sipation that is highly nonuniform, and distributed on a 
spatial multifractal. In contrast, dissipation range fluctu¬ 
ations are characterized by a linear ^ (m). This indicates 
global scale invariance and is associated with a distribu¬ 
tion of energy dissipation that is also nonuniform, but 
in this case distributed on a monofractal. Hence, the 
probability density functions (PDFs) of dissipation range 


magnetic field increments should collapse onto a unique 
scaling function. Figure 2(d) shows PDFs correspond¬ 
ing to T = {0.01,0.22,0.44,0.67,0.89,1.11,1.44} s that 
are rescaled by their standard deviations and overlaid, 
where the smallest r shows the associated errors. It is 
apparent that there is a very good collapse onto a single 
curve even up to several standard deviations. The largest 
events are not well-sampled as indicated by greater sta¬ 
tistical spread at larger values of the increments. A fitted 
Gaussian distribution illustrates the non-Gaussian PDF 
tails, which reflects the presence of rare large amplitude 
fluctuations. 

These results contradict earlier studies which found 
multifractal scaling at kinetic scales using PIG simu¬ 
lations [ 3 ^ and from direct observation of the scale- 
dependent kurtosis [s^. This inconsistency could be be¬ 
cause the simulations were set within the complex topol¬ 
ogy of a reconnection region whereas our observations 
sample the high speed outflow jet, and that kurtosis is 
sensitive to very large fluctuations which are not statis¬ 
tically well-sampled in heavy-tailed distributions. How¬ 
ever, the higher-order scaling shown in Fig. 2 is almost 
identical to that observed in the solar wind on both MHD 
and kinetic scales [13, Ej . This suggests that magnetic 
field fluctuations generated by reconnection exhibit a de¬ 
tailed correspondence with intermittent turbulence, in¬ 
cluding a cross-over from multifractal scaling to global 
scale invariance and distinct non-Gaussian statistics in 
the inertial and dissipation ranges. 

Turbulence cascades energy from structures on larger 
to smaller spatial scales. The corresponding spatial field 
is nonuniform with strong fluctuations that have non- 
Gaussian statistics. These fluctuations are small scale 
coherent structures which sup por t spatial gradients that 
can contribute to dissipation (39j |. In analogy with hy¬ 
drodynamic turbulence this suggests the possibility that 
reconnection, also related to activity on small scales, con¬ 
verts magnetic energy into kinetic and thermal energies. 
The connection between these dynamical processes is ex¬ 
amined within the earthward flow interval by computing 
Di = J • E, the work done by electromagnetic fields on 
the particles in the spacecraft frame. Although this is 
not strictly a measure of irreversible dissipation it must 
necessarily include the work done to convert stored mag¬ 
netic energy into heat. Indeed, the identification of Di 
as dissipation is complicated by contributions from par¬ 
ticle acceleration, fluid motion, field line stretching, and 
compressions. To avoid some of this ambiguity [e.g. jd^, 
the work done is evaluated in a frame moving with the 
bulk proton velocity Dp = J • (E -|- V x B). 

Figure 3(a) shows the broad and slightly asymmetric 
PDFs of the work done in the laboratory and proton 
frames rescaled by the standard deviation. These distri¬ 
butions are almost identical, which implies the contribu¬ 
tion to the work done from the convective electric field is 
minimal. There is an excess of positive values in the dis- 
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FIG. 3. PDFs of (a) spacecraft and proton frame work done, 
(b) current density and (c) electric field fluctuations. Gaus¬ 
sian fits (dashed curves) are applied. The (d) mean proton 
frame work done conditioned on local current density thresh¬ 
olds, where J/ Jrms is larger than and equal to some value n. 
Also, the fraction of data in these averages is plotted. 


tribution cores and the average work done is {Di) = 72±8 
pWm“^ and {Dp) = 62 ± 8 pWm“^. These may be 
interpreted as (imperfect) estimates of net magnetic en¬ 
ergy dissipation into plasma internal energy, but will also 
contain other effects. Note, values beyond the dashed 
vertical lines (±2.5 (t) are not statistically well-sampled 
and are likely dominated by unphysical fluctuations on 
scales smaller than the spacecraft separations. Nonethe¬ 
less, a fitted Gaussian distribution shows these heavy¬ 
tailed PDFs are manifestly non-Gaussian. This is a di¬ 
rect indication of intermittency of dissipation, which in 
the context of the KRSH, is associated with the presence 
of intermittent fluctuations of the fluid variables such as 
B. Figures 3(b)-(c) show the current density and electric 
held huctuations, which together constitute Dp, are also 
independently non-Gaussian and intermittent. 

If the work done in the proton frame is highly struc¬ 
tured as implied by its leptokurtic PDF, then this in¬ 
homogeneity should be evident in suitable conditional 
statistics. Figure 3(d) shows averages of Dp conditioned 
on thresholds of the local current density {Dp\J)/{Dp). 
The fraction of data points used in each of these condi¬ 
tioned averages is also plotted F{n) = Yl' f/Yl where 
Y' only includes points that satisfy the threshold condi¬ 
tion J/ Jrms > n. This represents a reasonable and eas¬ 
ily accessible measure of the volume hlling factor of the 
regions satisfying the corresponding condition. These di¬ 
agnostics illustrate explicitly the nonuniform and patchy 
character of the work done since the normalized condi¬ 
tional averages of Dp strongly increase with smaller vol¬ 
ume fraction. This results in a mean Dp for the threshold 
J > 2 A Jrms that is about 3-6 times the global average, 
despite such high local current density regions occupy¬ 


ing less than 2% of the data. Hence, regions of higher 
electric current density are increasingly rare, but make 
disproportionately large contributions to the total work 
done. 

Discussion .—We have used Cluster multispacecraft 
data to examine the structure within a previously iden¬ 
tified magnetic reconnection jet in the terrestrial magne¬ 
totail. In particular, the intermittent nature of the mag¬ 
netic field component increments 6Bj., SBy, SB^ and the 
work done, J-E, on particles by the electromagnetic fields 
is characterized. The significance of this result is seen by 
recalling the structure of the KRSH for hydrodynamic 
turbulence: Svi ^ . It is postulated that the lon¬ 

gitudinal velocity increments bvi on scale i are related 
statistically to the dissipation rate e/ averaged over a vol¬ 
ume of size (?. For a low density plasma such as the mag- 
netospheric plasma sheet, a formal statement of refined 
similarity has not been elucidated. However, our analysis 
is tantamount to a statistical examination of both princi¬ 
ple elements of a putative analogous relation for plasmas 
[e.g. 13, 4^, assuming the time increments employed here 


are comparable to spatial increments which is reasonable 
for a form of random sweeping. 

The magnetic field fluctuations (increments) within a 
reconnection jet exhibit a multifractal non-self-similar 
scaling of higher order moments in the inertial range, 
which transitions to a self-similar, but still non-Gaussian, 
monofractal scaling in the the kinetic range. In addi¬ 
tion, we find that J • E within the same jet is highly 
non-Gaussian, with heavy tails in the probability dis¬ 
tribution. Regions of strong J ■ E are non-space-hlling, 
with indications that the large transfer of random en¬ 
ergy to particles is almost certainly highly concentrated 
in small volumes that contain atypically large electric 
current density. Thus, even if the magnetotail plasma is 
not ohmic in nature, its dissipation is statistically asso¬ 
ciated with regions of high current density. This finding 
is consistent with results from the explicit examination 
of electromagnetic work in two- and three-dimensional 
collisionless plasma simulations [i , 3 ^ |: and from less 
direct inference in observations of the solar wind Q , mag¬ 
netopause 0 and magnetosheath [^. The conclusion 
seems increasingly certain that intermittent dissipation 
is as typical in large low density plasma systems as it is 
in high Reynolds number hydrodynamic turbulence. 

While our analysis focuses on turbulence within a 
magnetic reconnection outflow jet, the reconnection pro¬ 
cess itself occurs within a small diffusion region. The 
NASA Magnetospheric Multi-Scale mission will make 
high-resolution plasma and magnetic field measurements, 
and thus should allow the statistical nature of J • E to 
be probed within the proton and electron diffusion re¬ 
gions. This will improve our understanding of magnetic 
reconnection and turbulent dissipation. Indeed, as di¬ 
agnostics for examining intermittent dissipation become 
more well understood, we may anticipate a more refined 
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expectation may be emerging regarding the heating and 
intermittency that will likely be observed by the upcom¬ 
ing ESA Solar Orbiter, NASA Solar Probe Plus and pro¬ 
posed ESA THOR missions as they seek to understand 
how the solar corona is heated and the solar wind is ac¬ 
celerated, leading to the emerging structure of the entire 
plasma heliosphere. 
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